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We determine the modulational stability of standing waves with small group velocity in quasi- 
onedimensional systems slightly above the threshold of a supercritical Hopf bifurcation. The stability 
limits are given by two different long-wavelength destabilisation mechanisms and generically also 
a short-wavelength destabilisation. The Eckhaus parabola is shifted off-center and can be convex 
from below or above. For nonzero group velocity the Newell criterion, which near the cross-over 
from standing to traveling waves becomes a rather weak condition, does not determine the destabil- 
isation of all standing waves in one dimension. The cross-over to the non-local equations that are 
asymptotically valid near threshold is discussed in detail. Close to the transition from standing to 
traveling waves complex dynamics can arise due to the competition of counter-propagating waves 
and the wavenumber selection by sources. Our results yield necessary conditions for the stability of 
traveling rectangles in quasi-twodimensional systems with axial anisotropy and form a starting point 
for understanding the spatio-temporal chaos of traveling oblique rolls observed in electroconvection 
of nematic liquid crystals. 



I. INTRODUCTION 

In this work the coupled complex Ginzburg-Landau equations 

dtA + udyA = ^lA + {l + ib2)dlA - c\A\'^A - h\B\'^A, (1) 
dtB - udyB = fiB + {l + ih2)dlB - c\B\^B - h\A\^B, (2) 

describing the interaction of two counterpropagating waves arising via a supercritical Hopf bifurcation in one space 
dimension (see e.g. [0) are investigated. Here A and B are the complex wave amplitudes (or envelopes). The 
coefficients c and h are complex (we write c = Cr + ici etc.) with > and + ft-r > (see below), all other 
coefficients are real. The (dimensionless) control parameter /i (> above threshold) carries the system across the 
instability. Time and length are scaled such that the diffusion constant is 1 Uniform scaling of Eqs.(|l|,H) with 
^3/2 ^jjj |-,g obtained by introducing slow variables T = fj,t and Y = fi^^'^y, and scaling A, _B, and the group velocity u 
with fj,^/^. Strictly speaking this requires u to be small of order 0(//^/^). While we focus on this case, we also discuss 
the case u = 0{1) as well as the cross-over between the two cases. We expect that our results will give insight into 
the case u = 0{1) somewhat further above threshold. 

The spatially varying part of the physical fields u are given in terms of these variables as 

u = ^i/^e*'^'-* (^A{y, t)e''^P'yki{x, z) + B{y, t)e'^-P-y^i2ix, z)) + O(^) + c.c. (3) 

where uih is the frequency of the fastest-growing linear mode, which at threshold {fi = 0) is the Hopf frequency, Pc the 
critical wavenumber and £1^2(2;, z) eigenfunctions of the linear problem. The hatted quantities (e.g. t) are unsealed 
physical quantities. 

In the case hr > Cr the two waves suppress each other and one often needs to consider only a single equation 
describing a traveling wave train. In many cases the group velocity term can then be discarded by going into a 
moving frame. The remaining equation has been studied intensely. In particular it is well known that the plane-wave 
solutions with wave number po are stable with respect to long- wave side-band perturbations inside the band restricted 
by 

1 + — - 2-!^ ^ \ ' > 0, (4) 
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where Ftw = — p^/cr is the amphtude of the wave. Thus the band center solution po = is the last to lose 
stability when the Newell criterion Cr + b2Ci > becomes violated In addition, for strong dispersion short-wave 

side-band perturbations can be important |^,|^. This occurs only for wavenumbers away from the band center and 
when the Newell criterion is satisfied (c^ -I- 620^ > 0). It is not relevant in the following. 

The coupled equations become important when becomes smaller than (or only slightly larger than) and then 
the superposition of oppositely traveling waves, leading to standing-wave solutions, becomes important. The stability 
of standing waves with respect to side-band perturbations has been studied so far only in a few special cases [||0J|,^,B|. 
The case u = and hi/hr = Ci/cr has been investigated in quite some detail in the context of polarized lasers . 
There, however, the effects discussed below do not arise. For general values of the coefficients the solution in the 
center of the wave-number band has been investigated, for which the stability analysis is substantially simplified 
10,^. In addition to a long- wave instability an instability arising first at a finite perturbation wavenumbcr has been 
identified. The nonlinear evolution of the long- wave instability has been studied using coupled phase equations 
The evolution arising from the instability at finite wavenumber has been studied through numerical simulations of 
the Ginzburg-Landau equations [^. The behavior away from the band center has not been studied in detail. From 
previous work we know, however, that the instability of the band-center solution may not always refiect the behavior 
of the solutions with other wave numbers | [lO| | . We therefore present here a comprehensive analysis of standing waves 
for the complete band of wave numbers. 

We know of no physical (quasi-) onedimensional system in which standing waves have been clearly observed. 
In thermal convection of binary mixtures in porous media in the Hopf bifurcation range rather strong oscillations 
of the Nusselt number have been observed which can be taken as a hint for their occurrence in that system. 
Theoretically, a Hopf bifurcation to standing waves has been predicted for rotating convection at small Prandtl 
numbers p2[ |. 

In quasi-twodimensional systems with axial anisotropy an analogous situation arises when oblique rolls, i.e rolls 
appearing at an oblique angle with respect to the preferred axis, superpose to give traveling rectangles. In fact our 
investigation is motivated in particular by recent experiments on electroconvection in thin layers of nematic liquid 
crystals in the usual geometry with planarly aligned director, (see e.g. |ll3| ) where in some parameter range Dennin 
et al. 
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found near onset extended small-amplitude spatio-temporal chaos characterized by the interaction 
of obliquely traveling waves (rolls) and defects in these waves. Figure |^ shows a typical state It consists of 

patches of waves traveling obliquely to the direction singled out by the anisotropy of the liquid crystal and their 
superposition, which can be considered as traveling rectangles |]l7| . A starting point for the understanding of this 
state of spatio-temporal chaos is therefore the consideration of the stability of traveling rectangles. 



FIG. 1. Oblique waves 
nematic liquid crystals 111 



zigs 




and traveling rectangle patterns observed experimentally in electro-convection of 



The liquid crystal system is axially anisotropic and the complex dynamics arises immediately above the onset of 
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convection in a supercritical bifurcation. Combined with the fact that the origin of the travehng waves in the hquid- 
crystal system is finally understood (on the basis of a weak-electrolyte model these features promise 

that the chaotic state may be understood quantitatively using equations which to a large extent can be dealt with 
analytically, i.e. Ginzburg-Landau equations. 

A full, small-amplitude description of this system requires 4 complex modes for the left- and right-traveling waves 
in the 2 oblique directions. If one considers the common situation where oppositely traveling waves suppress each 
other strongly one is lead to two coupled Ginzburg-Landau equations, 

dtA + udyA = fiA + (1 + ibi)dlA + (1 + ib2)d^A + 2adlyA - c\A\'^A - h\B\^A, (5) 
dtB - udyB = fiB +{1 + ibi)dlB + (1 + ih2)dlB - 2adl B - c\B\^B - h\A\^B. (6) 



In a restriction to modulations in the y-direction (|1|,^ are recovered. In this paper we therefore capture a subspace 
of the possible instabilities of traveling rectangles. For a quantitative comparison with the experiments the full 
twodimensional equations using parameters determined from the hydrodyamic equations pl| ] will be needed [ p2[ . 
Presumably the onedimensional situation is also accessible experimentally by studying properly dimensioned narrow 
channels. 

In Sec. |l| we study the stability of standing- wave solutions of Eq. We find that, surprisingly, their stability 

is limited by two types of long-wave instabilities and a short-wave instability. One of the long-wave instabilities is 
related to the Eckhaus instability. Its parabolic stability limit is, however, shifted with respect to the neutral curve 
and can even be convex from above instead of below (see fig.^ below) . The shift results from the linear group velocity 
term. Moreover, the stable solutions can be outside the Eckhaus parabola instead of inside. Not too far from the point 
of degeneracy between standing and traveling waves the Newell-type criterion for the instability of the band-center 
can be satisfied even for weak dispersion, but it docs not indicate the loss of stability of standing waves at all wave 
numbers. 

Thus the group velocity u plays an important role. In the general case of m = 0(1) the asymptotic validity of 
Eqs.(0j2|) as — !■ has been questioned and in that limit equations have been proposed which are coupled through 
a term in which the counterpropagating wave is averaged over [ p3[|2^ , |25| , p^ . The long- wave stability properties of 
standing waves in these equations are found to differ markedly from those obtained from Eqs.(0,0). This has been 



demonstrated explicitly for the waves in the band center [g^. In Sec. Ill we show that the long-wave stability 
results of the non-local theory can be recovered from those of Eqs.(0,||) by taking the limit ^ before allowing the 
wavenumber p of modulations to become small [ p4| . This corresponds to the situation in a finite system. The usual 
long- wave stability theory of Eqs.(^|H), on the other hand, amounts to p ^ in the infinite system. Strikingly, the 
cross-over from the result of the non-local equations to that of can occur already very close to threshold. 



II. THE STABILITY OF STANDING WAVES 

The (unmodulated) standing- wave solutions of Eqs.(|l|,|^) are given by 

with 

^{pL- pl)/{cr + K), oj ^ -upo - 62P0 - (c» + h^)F^, (8) 

Clearly one needs c^. -\- h,. > in order that the standing waves bifurcate supercritically. The neutral curve is given 
by i^^ =0. Constant phases can be added to A and B separately. 

The stability analysis with respect to side-band perturbations is based on the expansion 

A= (l + vie^'+'Py + wae^**'*™) p^^{^t+poy)^ 

B=(l + wie^'+'Py + W2e^'*-'Py^ p^^i-^t-poy) _ (10) 

Inserting (p|,p^ into (|l|,^ yields after linearization in Vi and Wi a dispersion relation in the form of a quartic polynomial 
for A — X{p) for given wave number pq and control parameter /i. Alternatively, through (^), the amplitude F can 
be considered as the control parameter. First of all, for p = 0, one finds a doubly degenerated eigenvalue A = 
corresponding to the two phase modes, and two amplitude modes with A = —2{cr ± hr)F'^. Thus the conditions 



3 



Cr > and > h"^ are needed for amplitude stability. The long-wave expansion of the dispersion relation of the 
relevant phase mode 



p = €pi, A = eAi +e-^A2 + with e < 1, (11) 

can be treated analytically. We find non-trivial instabilities at orders 0(e) as well as O(e^). Our analysis of those 
instabilities shows that quite generally one has to consider also instabilities at finite values of p. This is in contrast 
to the stability of traveling waves {A — or B = 0) for which no instability occurs at order 0(e) and only in extreme 
cases a short- wavelength instability arises We have therefore also performed a general stability analysis by 

solving the quartic polynomial for A numerically. It shows that despite the importance of short-wave instabilities the 
long-wave analysis provides over a wide range of parameters a framework for the overall structure of the stability 
regions. 

In the long-wave expansion one finds at leading order in e 

^'-'P'cf^^ ^^^^ 

with 

ai,2 = u{Cr ± hr) + 2po{&2(Cr ± K) - (Ci ± hi)}. (13) 

Thus, to this order the solution in the band center, pq = 0, is always stable since Ai is purely imaginary. In previous 
work therefore no instability was found at this order. Such an instability arises, however, for po ^ (and u 0) 
when the factors in the numerator of Eq. (p^ ) change sign. As in a conservative system, the two imaginary eigenvalues 
then turn into two real eigenvalues, one positive and one negative, and lead to a steady instability. This occurs at 

„(1) ^ " Cr-K (1) ^ " Cr+K . . 

2 b2iCr - hr) + - C,' ^"'^ 2 62 (c. + /l.) " /l, " C, ' ^ ' 



Note, that to this order the stability limit is independent of F and therefore also of p.. From ( |l^ one finds that the 
standing waves are un 
and vice versa, where 



standing waves are unstable outside the interval ipl^l^p'^l) and stable (to this order) inside that interval for 71 < 



71 = (c, - 62C.)2 - {h, - b2hr)^ = "'4rr7?T- (15) 

Clearly this is equivalent to pil] and p^^l having opposite or same sign, respectively. The stability limits ( p^ are 
sketched by dotted lines in fig.@. 



(a) \ i 




/ / \ 




1^ 








Wave Number p / 

/ 






\ Wave Number p 



FIG. 2. Sketch of possible long-wave stability limits for standing waves. Standing waves are stable inside or outside the area 
delimited by the vertical dotted lines (instability at 0(e)), depending on whether they exclude or include the band center, and 
inside or outside the parabolic area given by the heavy line (0(e^)). The dashed line denotes the neutral curve. 
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At O(e^) one finds that either Ai = or 

— Ti2 J fl„ _ 

p2 '^■^ p2 



X2^pl\fi0-Pl^-^-P24l\ (16) 



with 



2b2hr{c,hr - Crhi) - Cricl - hi + - h"^) 

and given by (H). The limit po ^ {fi ^ 0) agrees with previous results 10,11 . The Eckhaus-stability limit A2 = 
is given by /i = ^'^'^'^^^ {po) with 

^(ecfe) = ^(c, + /i,)wpo + K'^'^'^'^pI where = 1 + ^^(c, + /i,). (18) 

Po Po 

In the parameter regime in which it is relevant (A^ < 0) the Eckhaus instability corresponds to an oscillatory instability 
and issues regarding its convective and absolute character arise. 
The stability limits given by (|6|) have three remarkable properties. 

i) At the band center pq —0 the standing waves become unstable for (3q > 0, i.e. for 

b2{hrh, - CrC,) + hi ~ cl > 0, (19) 

as already found in |Q. It was, however, not noted in that early work that this implies that near the transition from 
standing waves to traveling waves, hr ^ c~ , the solution at the band center becomes unstable for 

b2C^ < b2h„ (20) 

i.e. the condition involves only the imaginary parts 62, Ci, and hi, and even very weak dispersion can be sufficient to 
destabilize the standing waves in this regime. This is to be contrasted with the result for traveling waves for which the 
solution in the band center becomes unstable only if the dispersion is sufficiently strong, as expressed by the Newell 
criterion 62^ < — 1 ||,^. 

ii) The second striking feature is that the parabola giving the stability limit, /i — /x''^^'^^ (po), is shifted with respect 
to the neutral curve as sketched in fig.|[ This is due to the effect of the group velocity u. As a consequence, the 
instability of the solution at the band center does not imply the instability of all wave numbers (as one might have 
surmised in analogy with the case of traveling waves). Instead, when the solution with po — becomes unstable 
(cf. ([1^)), the apex of the parabola diverges to /z — > ±00 and the parabola merely flips over, becoming convex from 
above. Thus, an off-center band of wave numbers remains that is stable with respect to this instability (cf. fig.||b). 
Alternatively, when 

= (21) 

the parabola degenerates to a straight line and flips over in that way. In either case, for u 7^ 0, there is always a 
range of wavenumbers (in particular for either small positive or small negative po) over which the solution is stable 
with respect to the Eckhaus instability even on the neutral curve. In certain parameter regimes those solutions with 
vanishing amplitude are also stable with respect to the conservative instability at 0(e) and therefore stable with 
respect to all long-wave perturbations. By continuity, however, the solution cannot be completely stable since the 
trivial state F = is unstable. Therefore the long-wave analysis predicts that an instability has to arise which occurs 
first at a finite perturbation wave number. This is confirmed by the numerical evaluation of the dispersion relation 
X{p) (see below). 

iii) In contrast to the case of traveling waves and that of steady patterns near threshold the stable standing-wave 
solutions do not have to lie inside the parabolic Eckhaus stability limit; the case in which the stable solutions are 
outside the Eckhaus parabola occurs as well. Since A2 can change sign only on the Eckhaus parabola and the neutral 
curve, the condition that the waves are unstable outside the parabola is given by P2 < 0, i.e. 

2b2hric^hr - Crht) - Cr{cl - hi + cl - hi) < 0. (22) 
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In the absence of dispersion this is always the case (since Cr > hr). For sufficiently strong dispersion, however, A2 
can become negative for large po, and the waves are stable outside the parabola. Comparison with (|l5| ) shows that 
near the transition from standing waves to traveling waves, i.e. for hr c~ ^ the behavior of the two instabilities is 
complementary at large \po\: either the pattern is unstable due to Ai or due to A2. However, away from this transition 
both instabilities can arise for large po or neither. The latter case is particularly interesting since it implies that the 
pattern is stable with respect to long-wave perturbations even on the neutral curve. This parameter regime is limited 
by the conditions 71 = and /32 = 0. Solving for 62 and hi one obtains 

Acrit) _ Cr h^ — + h^ — , > 

"2 — h ~Z T. ^ ' 

and 

h^^''^^ ^{l/Cr){c,hr±{Cr~'hr)\c\) Or hf'^'^^{l/Cr){c^hr±{Cr + hr)\c\). (24) 

This leads to three intervals in hi over which A^ < and (32 — 0. For /i^ > they are given by (— oo, /i-''^), (/i.P\ 

and (/ip\oo). For /i^ < the values h^j^''^^ and h^^''^^ interchange their roles. In these intervals and for &2 close 

to b^^^^^ the long-wave stability analysis suggests that the pattern is stable even on the neutral curve. Again, this 
indicates the appearance of an instability at finite wavelength (see below). 
Thus, a number of qualitatively different situations can arise. For 

K^^'^''^ > (25) 

the parabola is convex from below and otherwise convex from above. In either case the solution in the band center is 
stable for 

l3o < 0. (26) 

If the parabola is convex from below this implies that the solution is stable inside the parabola, otherwise it is stable 
outside the parabola. In the former case it is worthwhile to consider the width of the Eckhaus-parabola. For 

Pof32 > (27) 

it is narrower than the neutral curve and wider otherwise. This implies that P2 determines the Eckhaus-stability of 
the solution with large \pq\ (cf. (^2])). In each of the cases the conservative instability at 0(e) can arise inside or 

outside the interval (Po^i,Po^2) determined by 71. 

In figs.^^ representative cases for the stability limits are shown as obtained from the long- wave results (|l^) and 
(|l6|), and by numerically solving the dispersion relation X{p). In fig.^ the parabola is convex from below and the 
Eckhaus-stable solutions are inside. For the parameters chosen, the conservative instability destabilizes the solutions 
in a vertical strip inside that parabola. When the parameters are changed so as to render the band center unstable 
the parabola becomes convex from above and one obtains the situation given in figj3|D (Note the enlarged wavenumber 
scale). Only the small area between the parabola and the vertical stability limit corresponds to stable solutions. Since 
the long-wave instabilities do not arise on the neutral curve for small negative po an instability at finite wavelength 
arises there. It is indicated by solid squares. 
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FIG. 3. Stability limits for standing waves. Long- wave expansion at 0(e) (dotted line) and O(e^) (solid line). The dashed 
line denotes the neutral curve. The symbols denote the stability limits as obtained from the numerical analysis of the dispersion 
relation A(p) with the solid squares marking an instability at finite perturbation wave number. The parameters are as follows: 
a) c = 1 - 0.4i, /i = 0.6 + 0.3i, 62 = 0.4, u = 0.1, b) c = 1 - 0.6i, ^ = 0.7 + 0.6i, 62 = 0.6, u = 0.3. 



In fig.^ the Eckhaus parabola is wider than the neutral curve. Since the 0(e)-instability destabilizes the solution 
only in the small range close to the band center the wavenumber band is limited by the finite- wavelength instability for 
large |po|- In fig-^b the stable solutions are outside the Eckhaus parabola and again perturbations at finite wavelength 
are important. 




-0.6 -0.4 -0.2 0.0 0.2 0.4 -0.5 -0.3 -0.1 0.1 0.3 0.5 

Reduced Wavenumber Reduced Wavenumber Pj 

FIG. 4. Stability limits for standing waves (cf. fig.^. The parameters are as follows: a) c = 1 -I- 0.3i, h — 0.5 -I- 0.5i, 
b'z = -1.96, u = 0.2, b) c = 1 + 0.3i, /i = 0.2 + 0.5i, 62 = -4.85, u = 0.2, 

For the parameters chosen in figj^a the stable solutions are also outside the parabola, but the parabola is now convex 
from above. In the final case, shown in fig.^jb, the standing waves are again Eckhaus-stable outside the parabola, 
which is convex from below (cf. inset). Now the conservative instability limits the wavenumber band on the other 
side (compare with fig.^) and only a very small region of stability remains. 




FIG. 5. Stability limits for standing waves (cf. fig.^. The parameters are as follows: a) c = 1 — 0.3i, h — 0.7 + 0.7i, 62 = — 1, 
u = 1, b) c = 1-1- 0.7i, h^0.9 + 1.5i, ^2 = 0.7, u = 0.5. 



As indicated above, due to the finite group velocity there is generically a range of wave numbers for which the 
solution is stable on the neutral curve with respect to long-wave perturbations. Figs.^^ confirm that therefore 
generically instabilities occur at finite perturbation wave number n. This is to be contrasted with the case of traveling 
waves where they appear only for sufficiently strong dispersion [y|5[| . 

Let us briefly discuss the case of vanishing group velocity separately. As is apparent from (O) and (111) the stability 
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with respect to long- wave perturbations does not depend on the wavenumber po to 0(e): either all wavenumbers are 
unstable (71 < 0) or stable (71 > 0). In the latter case the long-wave stability limits are determined by ([l^). At O(e^) 
one obtains a symmetric Eckhaus parabola. Again it can be convex down or up and can be wider or narrower than 
the neutral curve. Two cases are shown in fig.ga,b. In fig.^ the Eckhaus parabola is narrower than the neutral curve 
and the stability limit is indeed given by the long-wave instability. In fig.^D the parabola is wider and the stability is 
limited by a instability with finite p. 
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FIG. 6. Stability limits for standing waves with vanishing group velocity. Long-wave expansion at 0(e) (dotted line) and 
O(e^) (solid line), and instability at finite perturbation wave number (solid squares). The symbols denote the stability limits 
as obtained from the numerical analysis of the dispersion relation A(p). Neutral curve given by dashed line. Note that in c) all 
wavenumbers are stable according to the long-wave analysis, but they are in fact all unstable due to an instability at finite p. 
The parameters are as follows: a) c = 1 -f l.li, h = 0.6 + 0.5i, 62 = —0.5, it = 0, b) c = 1 + l.li, h = 0.6 -I- 0.5i, 62 = 9, it = 0, 

For u — the instability at finite p can be captured analytically near the band center because it arises there for 
p — O{po). To that end the perturbation wavenumber p and the wavenumber pq are taken to be of the same order, 
i.e. 

p = epi, pq = epoi, A = eAi + e^A2 + ... (28) 

At O(e^) one obtains Ai = 0. Note that (|l^) does therefore not apply and one has to go to higher order. At 0{e^) 
one finds 

,(1,2) 



X\'''> ^pil3o±piVWi (29) 



with 



^1 = (^2 _ ^2)2 Pi - 4 ( - 2 + - ^]Poi- (30) 
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The limit po — > of the previous result (|l^),(p^ for Ai_2 is regained by taking the limit pi — > in ( p^ , ^0| ). Thus, 
if the first term in (^9|) is negative and the prefactor of ppi (0) is positive y/Wi is imaginary in the long-wave 
limit and the waves with small po are stable. However, the first term in Wi becomes important when the modulation 
wavenumber becomes larger. Thus, for \pi\ ^ \poi\ > ( ^9|) yields an instability at a finite modulation wavenumber 
if f3o + -0) > 0, i.e. 

62 (ci — hi) + Cr — hr < 0, or b2{Ci + hi) + Cr + hr < 0. (31) 

Such a case is shown in fig.^ in which all wavenumbers are stable with respect to long-wave perturbations {(3o < 0, 
> 0, and 71 > 0), but all wavenumbers are unstable to perturbations at finite pi. Note, that the relevant 
perturbation wavenumber goes to in the band center. In the special case po — the stability conditions ( |3l| ) are 
therefore also obtained in a long- wave analysis H . 



III. ANALYSIS FOR LARGE GROUP VELOCITY 



Recently there has been some discussion of the fact that equations (0J|) are not obtained in that form when 
a direct expansion of the microscopic equations like the Navier-Stokes equations is performed near threshold 
P, ^3|j24|j25| , |26| , ^ , p8| . Since the group velocity of the waves is generically of 0(1) one is, strictly speaking, either 
led to consider a hyperbolic system without the second space derivatives in (|l|,||) 12^J2^ or to introduce two different 
time scales Ti = et and T2 = e^t |||,^,|^,|5 26 1. One then obtains already at quadratic order a solvability condition. 



dr^A + udyA^Q, (32) 
dr^B^udyB^Q. (33) 

Solving these equations in the frames co-moving with the respective waves {Y± =Y^ uTi) one then obtains at third 
order 

dr^A = + (1 + ib2)dl^A- c\A\'^A -h< \B\^ > A, (34) 
dr^B ^ fiB + {l + ib2)dY_B - c\B\^B -h < \A\^ > B, (35) 

with < ... > denoting spatial averages. These equations differ from (|l|,|^) in particular in their non-local cross-coupling 
terms between A and B: At the band center the linear stability with respect to long- wave perturbations of the 
standing waves has been investigated within these equations and a result quite different from that obtained in sec.|| 
has been obtained It has been pointed out that the difference between the two long- wave results is due to the 
fact that two limits have been taken: small amplitudes, i.e. — *■ 0, and long waves, i.e. p — > 0. These two limits do 
not commute. It has been argued that the asymptotically correct sequence is to take first /i — > 0, which leads to the 
non-local equations (|3^ , ^ ), and then p ^ 0. It is therefore of great interest to understand the connection between 
the two results and in particular the respective regimes of relevance. As long as no long-wave limit is taken the local 
(reconstituted) equations are expected to be valid. 

As discussed in |Q the stability result obtained from the non-local equations can be recovered from the local 
equations if the limit /x — is considered before taking the long- wave limit. To extend the analysis presented in |2j] 
away from the band center we consider the full 4*'*-order dispersion relation of perturbations of standing waves within 
the local equations taking the limit — > before considering the long- wave limit. Due to the parabolic character of 
the neutral curve the wave numbers have to be rescaled along with /i, 

M = ?/^M2, Po^VPoi, P^VPi, F^-qFi. (36) 
The growth rates are also rescaled, 

\ = iup + rj^K. (37) 

One then obtains at 0{rf) 

+ A {2pl - Aipoib2Pi + 2crF^) - (1 + bl){4pl^ - pl)pl 

-4i{Crb2 - Ci)piPQiF^ + 2{cr + c,b2)piFl = (38) 

Note that m, which is assumed to be 0(1), does not affect the (convective) destabilisation. Actually this dispersion 
relation coincides with that obtained for traveling waves, if -Fi is replaced by the traveling-wave amplitude. If one 
now expands in small perturbation wavenumbers pi , one obtains for the growth rate 
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A = 2ipoipi- 



-Ci + Crh 



Pi 



\c\' 



r3 F-' 



■Poi 



0{pI). 



(39) 



Thus, in this hmit no conservative instabihty arises and the stabihty boundaries are given by a paraboHc curve. Note 
that the term in brackets corresponds to the negative of the expression in Eq.(|4|). In fact, the destabihzation of aU 
periodic solutions is indicated by the loss of stability of the solution in the band center which occurs for the same 
parameter values (as for the traveling wave). 



&2Ci < 



(40) 



This is in strong contrast to the result found in (16) above, where for non- vanishing group velocity the Eckhaus curve 



was shifted with respect to the neutral curve. None of the complex behavior discussed above is obtained in 

Central to the difference between the two stability results is the long- wave limit. To understand the connection 
between them it is useful to consider a distinguished limit in which ji and the perturbation wave number p are both 
small of the same order. 



Pq 



P 



F = eFi 



A — iup + e A4 



(41) 



which differs from the previous limit (pq ) in the scaling of the modulation wavenumber p and that of the growth rate 
A. It turns out that the wavenumber po of the standing waves themselves has to be taken to be of O(e^) in order to 
obtain a single expression that reduces to the results of the local and of the non-local equations in suitable limits. At 
leading order (0(e^)) one obtains then 

A4 = pI {2Ffhr{Crhi - ahr) UP02 - (c^ - h^) [cl - hl+ b2{CrCi - hrhi)) F^ - u'^CriCr + b2Ci)pl} + 

+ iN^^P2 {2(c^ - hi) (62(0^ - hi) + hrhi - CrCi) F^pa2 + Up\ (2upo2Cr(Cr&2 - Ci) - b2hr{Crh, - Cihr)Fl)} (42) 

Af = {cl - hl)F^ + u^pIcI (43) 
For small amplitudes Fi (| 



A4 



reduces to 



-P2- 



2po2hr- 



uct 



Cihr rn2 , ■ I r, Crb2 
■Ff +ip2 \ 2po2 



uct 



Fn+o{Fr 



(44) 



This limit agrees to leading order with the result from the non-local theory (^9|). Note that due to the choice po — 0{e^) 
the wavenumber dependence of A in the non-local theory (cf. (|3^)) does not appear to this order. In the limit of large 
amplitudes one obtains from (E2h 



A4 



pI 



ip2 { 2p, 



/3o-/3: 



PQ2U 

b2{4- 



+ 



UP2 



2 hrb2{Crhi — Cihr) 



(c2 - /l2)2F2 



0(Ff4) 



(45) 



with Po and /3i defined in (^7|). Comparison with (|lj) and ( p^ ) shows that the limit ( ^l| ) also reproduces to leading 
order the result from the local theory. The imaginary part corresponds to an expansion of Ai from (|l^ ) for small po, 
noting that a term iup2 has been split off in . 

Thus, the distinguished limit (^) provides a connection between the two regimes and allows to give an expresssion 
for the parameters at which the cross-over occurs. Focussing on the band center, po = 0, the stability is determined 
by the second and third term in (^2|). A reasonable definition for the cross-over is therefore their ratio F, 



F = 



(cl - hl)(cl -hl + b2{CrCi - hrhi)) 



.(Cr -I- b2Ci) 



u2p2 



(46) 



The non-local theory is appropriate for F <C 1. For fixed values of the group velocity u, this corresponds to small 
amplitudes F'^ ^ fi. The cross-over to the regime in which the results of the local theory are recovered occurs near 
^c, which is defined via F = 1. Eq.(^) shows that depends not only on the group velocity u but also on the 
wavenumber p of the relevant perturbations and therefore on the size L — 2tt /pmin of the physical system under 
consideration, 



(47) 
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This dependence is illustrated in fig.^a. There the analytical results ( ]l8| , |39| ) are shown as well as the numerical result 
from the full dispersion relation for three system sizes corresponding to Pmin — 0.02, Pmin — 0.01, and Pmin — 0.005. 
The other parameters are c = 1 + 0.5i, h = 0.5 + 0.5i, u = 10, and 62 = 0. As expected, the non-local theory is correct 
for small /i. For larger n strong deviations arise for po > 0. They signify the cross-over to the result of the local 
equations. In agreement with (E^), in longer systems the cross-over occurs at smaller values of fi. This illustrates 
the fact that the difference between the two approaches is due to the fact that the limits /i ^ and p ^ are not 
interchangeable. In the limit p — > at fixed /i one obtains the result of the local equations, whereas for — > at fixed 
p the result of the non-local equations are obtained. It is noteworthy, that the cross-over occurs at quite small values 
of /Lt even though the group velocity is quite large {u = 10). At the same time the instability for po < is not captured 
at all within the long-wave limit of the local equations. As discussed in sec.|l[ within the long-wave limit the waves 
would be stable all the way to the neutral curve, implying that a short-wave instability arises. Interestingly, although 
this instability appears in the local equations as a short-wave instability, it is captured correctly in the long-wave 
analysis of the non-local equations. 
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FIG. 7. Comparison between stability results from non- local (l^jsHl) and local equations (|l|,|). Solid line from (|l|), dashed 
line from (|3^), dotted line gives the neutral curve. Symbols denote the stability limits from the full dispersion relation for 
Pmin = 0.02 (open circles), pmin — 0.01 (triangles), pmin = 0.005 (sohd circles), and pmin = 0.0005 (in c)). The parameters are 
c = 1 -f 0.5i, h — 0.5 + 0.5i, u — 10, and 62 = (in a)) and b2 = —1.9 (in b) and c)). 

Eq.(^6|) suggests that the cross-over to the result of the local equations occurs at a lower value of fi ii Cr + 0^62 
is small, i.e. close to the Newell criterion ( p0[ ) for the standing waves (within the non-local equations). This is 
demonstrated in fig.^. There the stability limits are shown for the same parameters as in fig.0a except for 62 which 
is here —1.9. Thus, + 0,^62 ~ 0.05. As in the case 62 = 0, the long-wave results of the local theory are adequate for 
long systems and larger fi (note the different scales as compared to fig.^) for pq > 0. The long-wave stability limits 
obtained from the non-local equations deviate, however, quite strongly even at very small values of fi. Strikingly, 
before the parabolic dependence of the non-local result takes over when fj. is decreased, additional deviations arise, 
which are due to the terms of 0{p'^) in (|39|). They arc important for /i < oc p^. Of course, if the system size is 
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increased these deviations become smaller, but at the same time the cross-over is shifted to lower values /ic, as well. 
More precisely, one has the scaling 



Hf (X IJc oc p for p ^ 1 . 



(48) 



Thus, for sufficiently large systems there always remains an intermediate regime in which the parabolic behavior 
from the non-local equations is relevant. It cauj- however, occur for values of fj, that are extremely small. This 
is demonstrated in fig.^c which shows a blow-uptl of fig.0b for Pmin = 0.0005. The usual local Ginzburg-Landau 
equations contain this result and capture in addition the rich possibilities that can arise for larger values of fi. 



FIG. 8. Comparison between stability results from non-local ( p4| , p5[ ) and local equations (^jy). Within the non-local equations 
all standing waves are unstable. Solid line from (^8|), dotted line gives the neutral curve. Symbols denote the stability limits 
from the full dispersion relation for pmin = 0.005 (solid circles). The parameters are c = 1 -|- 0.5i, h = 0.5 -|- 0.5i, u — 3, and 



A further interesting case is obtained if the Benjamin-Feir criterion Cr + &2Ci is satisfied within the non-local 
equations. In that case all standing waves are unstable immediately at onset. At slightly larger values of the control 
parameter, however, the waves can become stable. This is shown in fig.|[ Although with increasing values of the 
group velocity the gap of instability increases, it is seen that the restabilization can occur already for quite small 
values of the control parameter even for moderately large group velocity. 

Of course, at larger values of fi higher-order corrections will have to be included. Since the expansion leading to 
(|l|,|^) is only asymptotic (rather than convergent) firm predictions for finite values of /x cannot be made. However, if 
the cross-over to the local equations occurs for small values of /i it is quite reasonable to expect that the rich behavior 
displayed in figsj^||will not be strongly affected by the higher-order corrections. 



In order to determine the dynamics ensuing from the instabilities found analytically we have solved (y|^) numerically 
for a range of parameters. We have focussed in particular on the behavior in the vicinity of the transition from standing 
waves to traveling waves (hr — > c~), since in this regime even weak dispersion can lead to a destabilization of standing 
waves at all wavenumbers (cf. fig.0b) . Figures |9|-[ll| show a sequence of space-time diagrams for increasing values 
of hr straddling the transition point hr = Cr- As initial condition a slightly perturbed standing wave was chosen. 
The shades of gray in the space-time diagrams indicate the normalized difference {\A\ — |_B|)/(|A| + \B\) between the 
amplitudes of the right- and left-traveling wave with bright areas indicating domains in which the right-traveling wave 
dominates and dark ones those in which the left-traveling wave dominates. For hr = 0.7 both wave components are 
about equally strong resulting in a solution which corresponds mostly to a standing wave. This standing wave exhibits. 



'^For these parameter values round-off errors produced some scatter in the data. 




'0 



b2 = -2.2 (cf. fig.|). 



IV. NUMERICAL SIMULATIONS IN THE UNSTABLE REGIMES 
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however, complex dynamics due to the persistent occurrence of phase sHps in one of the two wave components and 
due to fairly locahzed disturbances propagating through the system. 




FIG. 9. Space-time diagram (space horizontally, time vertically upward) of evolution obtained by numerical simulation of 
(0,^) for /i = 0.4, c = 1 — 0.6i, h = 0.7 + 0.6j, 62 = 0.6, u = 0.3 (cf. figj^). The shades of gray indicate the normalized difference 
(|j4| — |_B|)/(jA| + |_B|) between the amplitudes of the right- and left-traveling wave (white corresponds to right-traveling, black 
to left-traveling waves). The length of the system is L = 400 and the time shown is from t — 4000 to t = 5000. 



Closer to the transition to traveling waves (/i^ — 0.9) the traveling-wave components become stronger and a 
characteristic behavior is seen (fig.|l^). Domains of left-traveling waves expand rapidly to the right for some time 
until the adjacent right-traveling wave starts to expand rapidly to the left pushing the left-traveling wave back. In 
the space-time diagram this sequence leads to dark triangles pointing to the right and bright ones to the left. Even 
for hr = 1.1, i.e. in the regime in which traveling waves of a given wavenumber are stable with respect to standing 
waves of the same wavenumber, these dynamics of rapid expansion and contraction persist as shown in fig.|l^. 
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FIG. 11. Space-time diagram (space horizontally, time vertically upward) of evolution obtained by numerical simulation of 
(0,^) for hr — 1.1. (other parameters as in fig.^. 

The origin of the complex dynamics shown in fig.|ll| is related to the selection of a non-zero wavenumber by the 
fronts in |^| and \B\, i.e. the sources, combined with the competition of counter-propagating waves of different 
wavenumbers j29j . A simple calculation shows that right-traveling waves with a wavenumber pg = pg'^' are stable with 
respect to left-traveling waves with wavenumber po = only for 



fir — I 



(49) 



Beyond this stability limit infinitesimal left-traveling waves are not suppressed by the right-traveling waves and grow. 
If the instability is only convective the left-traveling waves are swept away and, depending on the size of the domain 
of right-traveling waves, either perturb the next source or they grow to full amplitude forming a new domain of 
left-traveling waves. These dynamics ensue if the wavenumber selected by the source [ ^o| is beyond the stability 
limit given by (49). We have numerically determined the selected wavenumber. As shown in fig.^ it reaches the 
stability boundary ( ^9| ) at hr — 1.3. Thus, only for hr > 1.3 the wavenumber selected by a source is stable to all 
counter-propagating waves and the domains of left- and right-traveling waves become stationary and stable. These 
dynamics have also been discussed recently in the context of an investigation of sources |pl| . 
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FIG. 12. Wavenumber selected by source (squares) and stability limit of traveling waves with respect to counter-propagating 
waves (solid line). Between the dotted lines domains of traveling waves are unstable. (Parameters as in figJS). 
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V. CONCLUSION 



In conclusion, we have performed a detailed stability analysis of standing waves with small group velocity and found 
behavior that is quite different from that obtained for traveling waves. Due to the group velocity, the analogue of 
the Eckhaus parabola is shifted with respect to the neutral curve and can be convex from below and convex from 
above. The unstable solutions can be outside or inside the parabola. In addition, a conservative long-wave instability 
can arise already at lower order and generically a short-wave instability arises. The Newell- type criterion for the 
instability of the band-center can be satisfied even for weak dispersion, but it does not indicate the loss of stability of 
standing waves at all wave numbers. 

For group velocities of 0(1) the equations discussed here can be viewed as model equations that are obtained by 
reconstitution of the leading orders in the amplitude expansion. As has been shown previously [^,^ 25 2^ , in the 
asymptotic limit of small amplitudes non-local equations are obtained. Within the non-local equations the stability 
properties are considerably simpler. Our detailed analysis of the local (model) equations confirms that the non-local 
equations are valid immediately above threshold. However, in particular for large systems and for parameters close 
to the Newell criterion the control-parameter range over which the non-local equations are valid can be extremely 
small. Thus, already very close to threshold a cross-over can occur to the rich behavior discussed here. A very 
interesting question is therefore whether this behavior is observable in experimental systems, as might be expected 
if the cross-over occurs very close to threshold. In nematic liquid crystals the group velocities turn out to be quite 
small [^^0 . This sytem is therefore a good candidate to investigate this question. 

The stability analysis of standing waves gives sufficient conditions for the instability of traveling rectangles in 
two dimensions. It is in qualitative agreement with experimental results on spatio-temporally chaotic waves in 
nematic liquid crystals [p^, 15|jl^ ]. Based on the recent calculation of the nonlinear coefficients in the Ginzburg- 
Landau equations [ pTf , quantitative comparisons will be possible p^ . 

In view of the richness of the presented results experiments on rotating convection at small Prandtl numbers are of 
great interest, since in this system a Hopf bifurcation to standing waves has been predicted | p2[ . 

Finally we wish to emphasise that the phenomena treated in this work pertain to spatial modulations (in particular 
long-wave ones) of the two underlying traveling wave modes. The amplitude degeneracy occurring at hr = in 
principle leads to other effects and necessitates the inclusion of higher-order terms . 

Financial support by DFG (Kr690/12) and by DOE (DE-FG02-92ER14303) is gratefully acknowledged. 
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